Theory of point contact spectroscopy in electron-doped cuprates 
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In the hole-doped d x 2 _ y 2-wave cuprate superconductor, due to the midgap surface state (MSS), 
a zero bias conductance peak (ZBCP) is widely observed in [110] interface point contact spec- 
troscopy (PCS). However, ZBCP of this geometry is rarely observed in the electron-doped cuprates, 
even though their pairing symmetry is still likely the d x 2_ y -2-w&ve. We argue that this is due to 
the coexistence of antiferromagnetic (AF) and the superconducting (SC) orders. Generalizing the 
Blonder-Tinkham-Klapwijk (BTK) formula to include an AF coupling, it is shown explicitly that the 
MSS is destroyed by the AF order. The calculated PCS is in good agreement with the experiments. 

PACS numbers: 74.20.-z, 74.45. +c, 74.50.+r, 74.25.Ha, 74.25.Jb. 
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Pairing symmetry is an important issue towards the 
development of an understanding the mechanism of su- 
perconductivity. For hole-doped high-T c cuprate super- 
conductors, it is generally accepted that the pairing sym- 



metry is d x 2_ y 2-w&ve Among many supporting ex- 
periments, PCS measurement shows a ZBCP due to the 
existing MSS in the [110] direction @, S H- On the 
electron-doped side of cuprates, although no consensus 
has been reached yet, more and more recent experiments 
have found results also consistent with a d x 2_ y 2-w&ve 
pairing symmetry [j| [|| @, H, H, G3 ■ Thus one expected 
to see a similar ZBCP for the electron-doped cuprates. 
The situation is more subtle than this naive expectation 
however. The ZBCP has been consistently observed in 
the underdoped samples. In the optimally- and over- 
doped samples, in contrast, two dopin g-d ependent co- 
herence peaks were generally observed [111 . [l2l [T3L [l4| • 
This has been taken to indicate that the excitation gap in 
electron-doped cuprates might switch from d- to s-wave 
when the doping is increased - a phenomenon consis- 
tent with what was observed in the magnetic penetration 
depth measurement [l5j . 

The clues towards an understanding of the complicated 
PCS of the electron-doped cuprates lie in the doping evo- 
lution of the two (often called a and 0) pocket Fermi 
surfaces (FS) and the nonmonotonic d x 2_ y 2-w&ve exci- 
tation gap, as revealed by angle resolved photoemission 
(ARPES) measurements pi, [Tg|. These phenomena have 
been interpreted in terms of a phenomenological two- 
band model [13], which in turn led to a successful ac- 
count for the magnetic penetration depth measurement 
[ItJ and Raman scattering [l8[ data. The key feature 
of the two-band model is that a- and /3-band supercon- 
ducting (SC) gaps are both monotonic d 2 ,2_ l) 2-wave, but 
with a different (doping dependent) amplitude. Thus it 
was originally expected that [110] MSS should also exist 
in electron-doped cuprates. The difficulty in describing 
the experiments suggests that there may be some other 
physics intervening in the system in the SC state. 

The most promising candidate (for this absence) is the 
coexistence of an AF order with the SC order. This 
scenario finds support in ARPES as the data has been 
well explained in terms of a k-dependent band-folding 



effect due to an existing AF order [20, UM, \M \M- 
Apart from the ARPES data, no other direct evidence 
indicates that AF order does indeed exist in these ma- 
terials. Therefore, it is highly desired that some other 
experiments provide more definitive information regard- 
ing the possible existence of AF order in the electron- 
doped cuprates. When a normal electron is incident into 
a superconductor, it will induce the excitation of single 
quasi-particles (QPs) corresponding to all possible order- 
ing. The tunneling experiment is thus considered to be 
one of the best choices. Among various tunneling mea- 
surements, PCS is one of the most sensitive probes for 
the electronic states. 

In this paper, we shall focus on the effect of AF order 
on the PCS of a normal metal-insulator-superconductor 
(NIS) junction. It will be shown explicitly that the MSS 
can be destroyed by AF order. The observed PCS on 
electron-doped cuprates is actually the result of compe- 
tition between the contributions from AF and SC orders. 

Based on the simplest model that captures the essential 
physics, we will first consider a superconductor overlayer 
which is coated with a clean, size-quantized, normal- 
metal overlayer of thickness d, that is much shorter than 
the mean free path I of normal electrons. The interface is 
assumed to be perfectly flat and infinitely large. Consid- 
ering I — > 00, the discontinuity of all parameters at the 
interface can be neglected, except for the SC order pa- 
rameter to which the proximity effect is ignored Q. For 
d 3 .2_j,2-wave superconductor of interest, the interface is 
manipulated to be perpendicular to the k x axis (along 
the [110] direction). When both SC and AF orders ex- 
ist, QP excitations of an inhomogeneous superconductor 
have a coupled electron-hole character associated with 
the coupled k and k + Q [Q = (tt, 7t)] subspaces. Corre- 
spondingly QP states can be described by the generalized 
Bogoliubov-de Gennes (BdG) equations [24| 

Eu\ = houi + Ak^i + $^2 

Evi = AkUi — hoVi + $^2 

Eu 2 = $ui + h u 2 + A k+ QZ) 2 

Ev 2 = <JH>i + A k +QW 2 - h v 2 , (1) 

-ft 2 V 2 ./2r7i — fi with fi the chemical poten- 
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tial, <& is the AF order parameter, and Ak (= Ao sin 29) is 
the d x 2_ y 2-wave SC order parameter with Ak+Q = — Ak- 
Here the two-component wave functions u\ and v\ are re- 
lated to k subspaces, while u-i and V2 are related to k + Q 
subspaces. In arriving at Eq. ([1}, the pairing potential 
is assumed to be ~ Ak8 (x) with 9 (x) the Heaviside 
step function and Ak, given above, the Fourier trans- 
form of the Cooper pair order parameter in the relative 
coordinates. 

To better describe the effect of an AF order, the FS 
will be approximated by a square (see Fig. [T|). Thus at 
nearly half-filling, the FS matches the MBZ boundary. 
Under the WKBJ approximation pit] (/ = 1,2) 
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with 7 being the attenuation constant for |i? (kp) < 
A (kp)| and kp = (k x , k v , k z ), Eq. (P) becomes the An- 
dreev equation in the k x direction 
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for the superconducting overlayer (x > 0). Here e = 
e(fcx) = i'jkx/m. The wave-vector components parallel 
to the interface are conserved for all possible processes. 

Solving Eq. one obtains eigenvalues E = 

±yA^ + <f> 2 + e 2 where + (— ) corresponds to the 
electron- (hole-) like QP excitation. Since e(— k x0 ) — 
-e(k x0 ) and A (-k x0 ) = -A(fc x0 ), states for k x = k xQ 
and — k X Q are actually degenerate. Thus for k x — k x <j, one 
can have two degenerate eigenstates for electron-like QP 
excitation, while for k x — —k x o, one can have another 
two degenerate eigenstates for electron-like QP excita- 
tion. Superposition of these four eigenstates thus gives a 
formal wave function for the superconductor overlayer 
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Here E± = E ± e and a are the coefficients of corre- 
sponding waves. The wave function ips(%) remains the 
same if the eigenstates of hole-like QP are made up. 

Solving Eq. ([2]) when both Ak and $ are set to be 
zero, one may obtain bound states to the normal metal 
overlayer (— d < x < 0). In this case, the eigenvalues be- 
come E = ±k x oki/m assuming that the incident electron 
has the wave vector k x = k±. At the interface, the wave 
functions of normal metal and superconductor meet ideal 
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FIG. 1: Schematic plot shows all possible reflection and trans- 
mission processes for a normal electron incident into an NIS 
junction, where an AF order exists in the SC side. The inset 
shows the incident wave vector ki? = (fe^o, k y ,k z ) and its cor- 
responding AF wave vector kF + Q = k.FQ = { — k x o,k y ,k z ) 
due to an AF coupling. Both vectors are tied to the FS, which 
is approximated by a square (thick line). For convenience, k x 
axis is chosen to be along the [110] direction. 



continuity iPn(x = 0) = i^s{ x = 0). One thus obtains 
the formal wave function for the normal metal overlayer 



ip N (x) = 

/ e lklX A k \ 
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Considering the effect of free boundary at x = —d, it 
requires that iPn(% = —d) = 0. We then obtain the 
following eigencondition for the surface bound states: 



-2ik\ dj£ 



2ik\ d 



E_ = 2$. 



(5) 



Eq. ([5]) represents one of the major result in this paper. 
When the AF order $ = 0, there exists a zero-energy 
state which is responsible for the ZBCP widely observed 
in hole-doped d x 2_ y 2--w&ve cuprate superconductors 0]. 
When $ ^ 0, zero-energy state no longer exists such 
that the energy of the existing state is always finite. It 
is thus argued that the ZBCP does not exist in a system 
where the AF and SC orders coexist, and it is the origin 
which led to the absence of the ZBCP in the electron- 
doped cuprates. The above argument remains correct 
when d — > 0, which corresponds to the case of no normal 
metal overlayer (vacuum or with an insulating layer). 

In order to compare the theory with experiments, we 
next consider a superconductor [110] surface in point con- 
tact with a normal metal (STM) tip. In this case, a thin 
insulating layer is considered to exist between the normal 
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metal and the superconductor (see Fig. [I}. For this kind 
of NIS junction, the barrier potential is assumed to take a 
delta function, V(x) = HS (x). Considering that an elec- 
tron is injected into the interface from the normal metal 
side (with an angle 9), four possible reflections and their 
related coefficients are detailed as follows: (a) Normal 
reflection (reflected as electrons) with the coefficient rs- 
(b) Andreev reflection (reflected as holes, due to electron 
and hole coupling in the k subspace) with the coefficient 
Ta- (c) AF reflection (reflected as electrons, due to the k 
and k + Q subspace coupling) with the coefficient r e AF . 
(d) Andreev- AF reflection (reflected as holes, due to elec- 
tron and hole coupling in the k + Q subspace) with the 
coefficient r\ F . In terms of these coefficients, the for- 
mal wave function for the normal-metal side can then be 
written as 



exp (ik x0 x) + r s exp (-ik x0 x) 
r A exp (ik xQ x) 
r e AF exp (ik x0 x) 
r\ F eyip(-ik x0 x) 



The four reflection coefficients are determined by the 
boundary conditions: 



tpN (x) \ x=0 - = 1p S (x) \ x=0 +, 

2mH dips [x) 

-ips (x) | x =o+ = — — — U=o+ 



(6) 
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The normalized tunneling conductance is then given by 



a(E,9) = l-\r s (E,9)\ 2 + \r A {E,9)\ 2 



+ Kf {E,l 
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Actual tunneling conductance is intimately determined 
by the junction properties. When the tip and the su- 
perconductor are in an ideal point contact, i.e., when 
the effective potential barrier Z = 2mH/h 2 is sufficiently 
low and narrow, the wave functions at the interface meet 
good condition of continuity. Consequently Andreev and 
AF reflections are important. For non- ideal point con- 
tacts, in contrast, only the normal reflection is impor- 
tant. In the following, we focus on the case of low and 
narrow barrier height, i.e., the PCS. 

Since the theory is aimed to the electron-doped 
cuprates, the issues concerning their doping dependent 
AF and SC orders are crucial. Shown in the inset (I) of 
Fig. (Ha) is a typical SC gap for electron-doped cuprates. 
The gap is piecewise: both segments are fitted to the 
monotonic d x 2_ y 2-wa,ve (Aq sin 29), with different ampli- 
tudes (A 0/ 3 and A 0q ). The (3- and a-band FS segments 
are characterized by the cutoff angle 9\ and 62- The val- 
ues of doping-dependent 6% and 9i were extracted from 
the ARPES data 0. 

In a real NIS junction experiment, the total tunnel- 
ing conductance is given by o(E) — J o(E,9)d9. Here 
the integration over the angle between 9\ and 92 is ruled 
out, to which the FS is absent. In the following the di- 
mensionless quantities 3>g = $/A 0) 9 and <& Q = $/A 0a , 



TABLE I: Three (dimensionless) fitting parameters: <J?/3, $ a , 
and Z, obtained in Fig. [2] A0/3, Ao a , and $ (in unit of meV) 
in the right three columns are determined exclusively from 
$/3 and $ a together with the actual peak energy. 
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FIG. 2: (Color online) Fitting to the low-resistance G-V data 
of electron-doped cuprates. Data in frames (a),(c),(d) are 
taken from Ref. on Pr 2 - 2 ,Ce a; Cu04 at T = 1.43K, while 
data in frame (b) is taken from Ref. p3] on Nd2- a; Ce a ,Cu04. 
Fitting parameters are summarized in TABLE [I] Inset (I) 
in (a) shows a typical piecewise d x 2_ y 2 -wave gap, with 9i 
(82) characterizing the cutoff angle for f3 (a) band FS. Inset 
(II) shows the Z — > calculation, emphasizing the zero-bias 
"peak" structure in (a). 



which play the determining role on the actual PCS, will 
be given from the fitting calculations. The actual values 
of Aoq, and Aop can also be obtained through the 
fitting processes. Moreover in the calculation, the QP 
energy E will be replaced by E — iT with T characteriz- 
ing the finite lifetime of the QPs. The parameters from 
the best fittings are summarized in TABLE Q] As a mat- 
ter of fact, only parameters associated with the /3 band 
(i.e., associated with the nodal region) are sensitive to 
the fitting. 

Figure [Ha) compares the theoretical calculations with 
the conductance-voltage (G-V) curve of PCS data [26[ 
on underdoped Pr 2 _ :c Ce 2; Cu04 (PCCO) (x = 0.13). For 
this case, Z = 0.05 is used to simulate a low resistance 
(R = 9.8ft). The result is shown to be in good quan- 
titative agreement with the data, to which (&p,$ a ) = 
(0.1,0.2) are obtained. The smallness of Q a and in- 
dicates that the effect of AF order is inappreciable. As a 
consequence, a strong ZBCP appears, revealing the dom- 
inance of the d x 2_ y 2-wa,ve symmetry of the gap. This is 
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in full support to the conclusion drawn in Ref. [2(| that 
the paring symmetry for the underdoped samples is con- 
sistent with d x 2_ y 2-wa,ve. To explore the AF effect in 
more details, the G-V curve is recalculated for an ideal 
point contact (Z, T — ► 0). In this limit [see inset (II) in 
Fig. Ufa)], the zero-bias "peak" is actually a mix with 
a narrow plateau. The plateau, that resembles the fea- 
ture of an s-wave superconductor, arises simply because 
a small but finite AF order exists near the nodal region. 

The evolution of the AF order can be understood from 
the doping-dependent FS. In underdoped electron-doped 
cuprates, due to their more distant FS from the MBZ, the 
scattering about Q = (jr, tt) is weak. This results a small 
AF order. With increasing the doping, FS approaching 
and crossing to the MBZ, the scattering about Q becomes 
more important which leads to a more important AF 
order |2ll . [27j . Upon further increase of the doping, long- 
range order is destroyed such that the AF order decreases 
and vanishes eventually. 

A similar plateau (with a dip) structure has also been 
observed in PCS on underdoped Nd2- :c Ce :r Cu04 [Hj]. To 
fit this set of data, ($g, $ Q ) = (0.7, 0.9) and Z = 0.25 are 
obtained [see Fig.[2^b)]. Comparing the result with that 
in inset (II) of Fig. [D(a), the G-V curve changes from a 
plateau to a two-peak structure. The latter results from 
much higher 3>g and <I> Q and a slightly higher resistance 
Z. The dip at zero bias is an evidence that the ZBCP 
does not exist in this case. 

Consider more closely how the FS segments evolve as 
the doping changes. At low doping, a-band FS first 
emerges in the antinodal direction before the supercon- 
ductivity sets in. When the doping is increased, /3-band 
FS appears simultaneously with the appearance of su- 
perconductivity. Since MSS is the signature for a su- 
perconductor which has symmetrically a positive and a 
negative portion of the gap, as long as $ is inappreciable, 
the ZBCP can still be observed in underdoped samples 
no matter how the FS segments emerge, or even in the 
(unrealistic) case without the nodal (/3-band) FS. 



PCS with high resistance Z has been observed on 
optimally-doped PCCO (x = 0.15) [H]. In the best fit- 
ting [Fig. Hfc)], large ratios ($0, $ a ) = (0.75,1.0) and 
Z = 1.5 (for R = 18f2) are obtained. The strong effect of 
the AF order $ leads to a clear two-peak feature, consis- 
tent with the case of a dominant s-wave gap. Neverthe- 
less, theoretical curve deviating from the data at higher 
biases confirms that a higher Z is in use and the wave 
functions at the interface do not meet good condition of 
continuity. 

In Fig. [2jd), fairly good fitting is also made with the 
PCS data on overdoped PCCO (a; = 0.17) 0, to which 
large (&f3,$ a ) = (0.7,0.75) and a relatively smaller Z 
(for R = 2.60) are used. Again, a two-peak feature con- 
sistent with an s-wave gap manifests the largeness of <fr. 
In Ref. the G-V curve in Fig.[2fd) was also fitted us- 
ing the BTK model. Various pairing models were tested 
and it is the <i -Ms-wave symmetry that leads to the best 
result. This does support the scenario discussed in the 
current paper that single particle excitation is gapped by 
both SC and AF orders, A off = |A d + i$ = y/A 2 d \ + $ 2 . 
The actual PCS is indeed the competitive result between 
the SC and AF contributions. 

In summary, PCS of the electron-doped cuprates is 
investigated. It is shown explicitly that MSS of the d- 
wave superconductor can be destroyed by the presence 
of an AF order ($). Due to the smallness of $, ZBCP 
occurs in the underdoped sample, consistent with the 
d a .2_ y 2-wave pairing. A more important effect of $ results 
a two-peak feature in optimally and overdoped samples. 
The phenomenon that the peak energy first increases and 
then decreases as doping increases is soundly explained. 
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